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ABSTRACT 

We report on the first measurements of the isotopic ratio 14 N/ 15 N in N 2 H + 
toward a statistically significant sample of high-mass star forming cores. The 
sources belong to the three main evolutionary categories of the high-mass star 
formation process: high-mass starless cores, high-mass protostcllar objects, and 
ultracompact Hll regions. Simultaneous measurements of the 14 N/ 15 N ratio in 
CN have been made. The 14 N/ 15 N ratios derived from N 2 H + shows a large spread 
(from ~ 180 up to ~ 1300), while those derived from CN are in between the value 
measured in the terrestrial atmosphere (~ 270) and that of the proto-Solar nebula 
(~ 440) for the large majority of the sources within the errors. However, this 
different spread might be due to the fact that the sources detected in the N 2 H + 
isotopologues are more than those detected in the CN ones. The 14 N/ 15 N ratio 
does not change significantly with the source evolutionary stage, which indicates 
that time seems to be irrelevant for the fractionation of nitrogen. We also find 
a possible anticorrelation between the 14 N/ 15 N (as derived from N 2 H + ) and the 
H/D isotopic ratios. This suggests that 15 N enrichment could not be linked to 
the parameters that cause D enrichment, in agreement with the prediction by 
recent chemical models. These models, however, are not able to reproduce the 
observed large spread in 14 N/ 15 N, pointing out that some important routes of 
nitrogen fractionation could be still missing in the models. 
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Introduction 


Nitrogen is the fifth most abundant chemical element in the universe, and exists in the 
form of two stable isotopes: 14 N and 15 N. In the terrestrial atmosphere (TA), the typical 
atomic composition 14 N/ 15 N as derived from N 2 is ~ 270. This value is quite consistent 
(within a factor 2) with that measured in cometary nitrile-bearing molecules (e.g. Bockelee- 


Morvan et al. 12008 . Manfroid et al. 20091 ) and in primitive chondrites (e.g. Briani et ah 2009 


Bonal et ah 20101 ). but it is significantly smaller than the value measured in the present- 
day Solar wind, which is_^_441 ± 6 , as derived from the particles collected by the Genesis 


spacecraft (Ma rty et ah 20111). This value is comparable to that of the Jupiter’s atmosphere 
(Fouchet et ah 120041) . and it is considered to be representative of the so-called ’’proto-Solar 
nebula” (PSN), from which our Sun has formed. However, to understand how these values 
are linked to the initial chemical composition of the molecular environment in which the 
System itself was formed, and put constraints on chemical models, measurements of the 
14 N/ 15 N ratio in dense molecular star-forming cores similar to that in which the Sun was 
likely formed are mandatory. 

In dense molecular cores, where the direct observation of N 2 is impossible, the 14 N/ 15 N 
isotopic ratio is measured from abundant N-bearing species such as N 2 H + , NH 3 , CN, HCN, 
and HNC, and the results are still puzzling. In fact, while the 14 N/ 15 N ratio measured 
toward pre-stellar cores or protostcllar envelopes in hydrogenated nitrogen compounds (such 
as N 2 H + and NH 3 ) is comparable to the value found in the PSN (e.g. in NH 3 ; Lis et 
al. 20 
et al. 


Gerin et al. 12009 : Daniel et al. 2013j), or even larger (1000±200 in N 2 H + ; Bizzocchi 
2013), the 14 N/ 15 N ratio derived from the nitrile-bearing species HCN and^HNC is 
significantly lower (140- 360, Hily-Blant et al. 2013a : 120-400, Adande & Ziurys 2012 : 160- 


290, Wampfler et al. 2 0141) . Differences in nitrogen hydrides and nitrile-bearing molecules 
within clouds are predicted by chemical models inclusive of spin-state chemistry (Wirstrom 
et al. 2012; Hily-Blant et al. 2013af) . but no models have been able so far to reproduce the 
low 15 N fractions measured in N 2 H + . The 15 N fraction measured in CN toward starless cores 
has been found marginally consistent with the proto-Solar value of 441 (500±75, Hily-Blant 
et al. l2013bh . in contr ast w ith the large 15 N fractions measured in HCN in both starless 
cores (Hily-Blant et al. l2013bl ) and protostcllar objects (Wampfler et al. 2014). This is also 


expected from theory (Hily-Blant et al. 2013bi). A gr adient of 14 N/ 15 N with the Galactic 
distance has been measured by Adande & Ziurys ( 2012 ) via observations of molecular clouds 
in CN and HCN, and found to be in agreement with predictions of Galactic chemical evolution 


1 Based on observations carried out with the IRAM-30m Telescope. IRAM is supported by INSU/CNRS 
(France), MPG (Germany) and IGN (Spain). 
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models (e.g. Clayton 2003], Romano & Matteucci 2003]). However, the 15 N fractionation of 
CN and HCN are often based on observations of the optically thin isotopologues containing 
13 C, because the main isotopologues are typically optically thick, which makes it less accurate 
the derivation of the column densities (e.g. Hily-Blant 2013a, 2013b). Roueff et at. ( 20151 ) 
have found possible reduced abundances of the 13 C due to the fact that nitriles and isonitriles 
are predicted to be significantly depleted in 13 C (Roueff et at. 120151) . However, this depletion 
is estimated to be of the order of the ~ 20%, namely often of the order of the calibration 
uncertainties associated to the measured column densities. 

In star forming regions, D- and 15 N- fractions are not correlated. For example, in 
L1544, NoH + is highly deuterated (D-fraction ~ 30%; Caselli et al. 2002), while an anti¬ 


fractionation is observed in 15 N (Bizzocchi et ah 120131) . Wirstrom et al. (2012) reproduced 
these differentiations and noticed that they are also present in meteoritic material, as only 
in a few cases 14 N/ 15 N correlates with H/D ratios (Mandt et al. 2014). This suggests an 
interstellar heritage of the Solar System 15 N and D isotopic anomalies. Therefore, it is 
important to gather more data and to put stringent constraints on current chemical models, 
with the ultimate goal of investigating the possible link between cometary and interstellar 
medium material (see also Caselli & Ceccarelli 2012), especially in intermediate- and high- 
mass dense cores^ since growing evidence is suggesting that our Sun was born in a rich stellar 
cluster (Adams 2Qlol) . 


In this letter we report the first observational measurements of the 14 N/ 15 N ratio in a 
sample of dense cores associated with different stages of the high-mass star formation process, 
and already studied in deuterated molecules. Our study can investigate if the nitrogen frac¬ 
tionation process varies with the physical evolution of the core. What is more, the 14 N/ 15 N 
ratio is derived from two well-studied nitrogen-bearing species, N 2 H + and CN, allowing us 
to investigate whether the fractionation of the two species follows different pathways as 
suggested by Hily-Blant et al. fl2013al) . 


2. Observations 


The spectra presented in this work are part of a spectral survey performed with the 
IRAM-30m Telescope in two millimeter bands, one centred at about 3 mm covering fre¬ 
quencies in the range 89.11 - 96.89 GHz, and another centred at about 1.3 mm covering 
frequencies in the range 216.0 - 223.78 GHz, towards 26 dense cores divided into three main 
evolutionary categories: High-mass starless cores (HMSCs); high-mass protostcllar objects 
(HMPOs); Ultracompact Hu regions (UC Hlls , see also Fontani et al. 201 ll for details on 
the source selection criteria). In Fontani et al. 020151) we have published part of these ob- 
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servations, namely that about the lines of CH3OH and of its deuterated forms, so that we 
refer to Sect. 2 and Table 2 of that work for a description of the observational parameters 
and technical details. In the 3 mm band, the fundamental rotational transition of N 2 H + , 
15 NNH + and N 15 NH + (at 93173.4 MHz, 90263.8 MHz and 91205.7 MHz, respectively) were 
observed simultaneously to the CH 3 OH observations described by Fontani et al. (120151 ) with 
velocity resolution of ~ 0.65 km s _1 (the telescope HPBW is ~ 26"). In the 1 mm band, 
several hyperhne components of the transition N — 2 — 1 of radicals 13 CN (50 components 
with rest frequencies in between 217.0326 and 218.0311 GHz) and C 15 N (10 components with 
rest frequencies in between 219.38749 and 219.93482 GHz) were also observed with velocity 
resolution of ~ 0.26 km s _1 (the telescope HPBW is ~ 11"). 


3. Results 


3.1. 15 NNH+ , N 15 NH+ , and N 2 H+ 


The 15 NNH + (1-0) has been detected in 14 cores: four HMSCs, six HMPOs and four 
UC HIIs. Because the separation between the hyperhne components is comparable to the 
line widths, the hyperhne structure is unresolved in almost all spectra. Therefore, the lines 
have been ht with single Gaussians, except in four cases (05358-nmi3, 05358-mml, ONI, 
19410) in which the ht to the hyperhne structure has given reliable results despite the small 
separation of the components and the poor spectral resolution. This simplified approach 
could affect the determination of the line width, but will not affect the measurement of 
the column densities because they will be derived from the total integrated area, as we will 
describe below. The N 15 NH + (1-0) has been detected in 11 cores: two HMSCs, hve HMPOs 
and four UC HIIs. The hyperhne components in this case are only three, with a faint central 
one and two fainter satellites, but the faintness of the satellites has forced us to ht many of 
the lines with Gaussians in this case too. Finally, the N 2 H + (l-0) line is clearly detected in 
all sources and the ht to the hyperhne structure gives good results in all spectra, pointing 
out that no hyperhne anomalies are present. In Fig. |T] we show the spectra of the 15 NNH + , 
N 15 NH + and N 2 H + (1-0) lines in the representative core 05358-nmi3. 


The total column density of N 2 H + , averaged within the beam, has been evaluated from 
Eq. (Al) in Caselli et al. (120021) . The formula assumes the same T ex for all the hyperhne 
components. For sources with total opacity well-constrained, T ex was computed from the 
output parameters obtained from the hyperhne ht procedure (in particular the total optical 
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depth, r t , and T ant x r t ) following the approach described in the CLASS user manual^. For 
the others, we have assumed the average T ex derived from the lines with well-constrained 
opacity in each evolutionary stage (18, 43 and 57 K for HMSCs, HMPOs, and UC Hlls, 
respectively), and obtained the column density from Eq. (A4) of Caselli et al. (120021 ). valid 
for optically thin lines. The spectral parameters used in our calculations (S'/i 2 , Einstein 
coefficients, rotational constants, energies and statistical weight of the_ levels) were taken 
from the Cologne Database for Molecular Spectroscopy (Muller et al. 20051 ). 


The total column densities of both 15 NNH + and N 15 NH + have been computed in all 
detected sources from the total line integrated intensity with Eq. (A4) of Caselli et al. (120021 ). 
assuming optically thin lines and LTE conditions. As excitation temperature, we have used 
that derived from the N 2 H + (l-0) line. All column densities and the parameters used to 
derive them (line integrated intensities and T ex ), are given in Table [0 


3.2. 


13 


CN and C 15 N 


Several hyperfine components of the iV = 2 — 1 transitions of both 13 CN and C 15 N 
were detected towards 13 cores, which yields a detection rate of 50%. As expected, most of 
the cores were detected in the hyperfine components with larger values of S/j, 2 , which are: 
J,F 1 ,F= (5/2, 3,4) - (3/2, 2, 3); (5/2, 3, 3) - (3/2, 2, 2); (5/2, 3, 2) - (3/2, 2,1) for 13 CN , 
and J,F— (5/2,3) — (3/2,2); (5/2,2) — (3/2,1) for C 15 N. Representative spectra obtained 
towards core 05358-nnn3 centred on these hyperfine components are shown in Fig. [0 

To compute the column densities of both radicals, we have derived first the integrated 
intensity of the group of hyperfine components (Table [ 2 ]) mentioned above, and computed 
the total column densities from Eq.(A4) of Caselli et al. (120021 ). which assumes the same 
excitation temperature in the various rotational levels, and it is valid for optically thin lines. 
The formula has been adapted to take the hyperfine structure into account, namely we have 
adopted the spectroscopic parameters (Einstein coefficients, energies and statistical weight 
of the levels) of the hyperfine components considered only. All spectroscopic parameters, as 
well as the partition functions used to derive the total column density, have been taken from 
the Cologne Database for Molecular Spectroscopy (Muller et al. 120051) . When the hyperfine 
components are not blended, their relative intensities are consistent with the assumption of 
optically thin lines and no hyperfine anomalies (within the calibration uncertainties). As 
excitation temperature, we have used the T ex reported in Tabled derived from N 2 H + . This 
assumption is critical, but the C 15 N/ 13 CN column density ratio does not change by varying 


2 for details: http://iram.fr/IRAMFR/GILDAS/doc/html/class-html/class.html/ 
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T ex of even an order of magnitud e, as it was already noted in other works dealing with 


isotopic ratios (see e.g. Sakai et al. 2012|, Fontani et al. 2014aJ). 


The CN column densities were then computed from the 13 CN ones by applying the 
isotopic ratio 12 C/ 13 C derived from the relation betw een this ratio and the source galacto- 
centric distance found for CN by Milam et ah 0200511 . Galactocentric distances have been 
taken from Fontani et ah ( 2014bl ). We have not corrected the column densities for a filling 
factor because we do not have any direct estimate of the size of the emitting region of both 
molecules, and we do not see any physical or chemical reason by which the emission of 13 CN 
and C 15 N should arise from different regions. 


In 19410+2336, and ONI, the 13 CN hyperfine components mentioned above are barely 
detected, while in G5.89-0.39 they are heavily affected by a nearby CH 3 OH stronger line. 
Therefore, in these sources we have used the better detected hyperfine component J,Fi,F = 
(3/2,1,1) — (1/2,0,1). In G5.89-0.39 the blending problem affects also the C 15 N lines, so 
that we have used the isolated hyperfine component J,F— (3/2, 2) — (1/2,1) to compute the 
C 15 N column density. For I00117-MM1, 121307, I19035-VLA1, and 23033+5851, detected 
in 13 CN but not in C 15 N, we have derived an upper limit for the column density of C 15 N. 
All column densities, and the line integrated intensities used to compute them, are listed in 
Table [2] 


4. Discussion 

The comparison between the column densities of the 15 N-containing species and those of 
their main isotopologues, derived as explained in Sects. 13. II and 1X21 is shown in Fig. [21 The 
corresponding 14 N/ 15 N isotopic ratios are given in Tables Q] and [2j The first result emerging 
from these numbers and from Fig. [2] is that the NNH + / 15 NNH + and NNH + /N 15 NH + ratios 
vary between 180 and 1300, namely over a range much larger than observed so far towards 
low mass objects in any N-bearing species. Some values are, therefore, compatible with both 
those of the TA and PSN, but much larger values are also measured. In addition, there 
is no evidence of a trend of these ratios with the evolution, because Fig. [2] shows that the 
sources belonging to different evolutionary stages possess similar 14 N/ 15 N ratios, within the 
uncertainties. Consistently, we do not find any correlation with the line widths, which are 
thought to be good indicators of the evolutionary stage. Therefore, whatever is the reason 
for this different 15 N enrichment, time does not seem to play a role. 

The second result is that the CN/C 15 N ratio is within ~ 190 and ~ 450, namely con¬ 
sistent within the errors with the 15 N enrichment of both TA and PSN. Interestingly, when 
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Table 1: Integrated intensity of the 15 NNH + , N 15 NH + and N 2 H + (1-0) transitions, total 
beam-averaged column densities computed as explained in Sect. 13.11 and corresponding 
14 N/ 15 N isotopic ratios. Uncertainties in the isotope ratios have been computed from the 
propagation of errors on the column densities. These latters have been derived from the 
errors on the line areas for the optically thin lines, given by rms x At> x s/N (rms = root 
mean square noise in the spectrum, Av = spectral velocity resolution, N = number of chan¬ 
nels with signal). For the N 2 H + (1-0) lines with well-constrained opacity, w e prop agated the 
errors on Tt and T ant x r t . The source names are taken from Fontani et al. 120111. 


source 

15 nnh 

+ (1-0) 

N 15 NH+ (1-0) 

N 2 H+(l-0) 

n 2 h+ 

n 2 h+ 

^ex 


/ TmbJ" 

!V( 15 NNH+ ) 

J i mbJ» 

1V(N 15 NH+ ) 

/ iMBdt 1 

JV(N 2 H+) 





K km a -1 

(Xl0 10 cm- 2 ) 

K km s 1 

(Xl0 10 cm- 2 ) 

K km s 1 

(Xl0 13 cm“ 2 ) 



(K) 

HMSC 

100117—MM2 

0.029(0.001)'' 

5.68(0.09) 

0.04 

< 6.8 

3.966(0.007)“ 

3.8(0.5) 

670(98) 

> 561 

19 

AFGL5142-EC™ 

0.040(0.002) 

24(1) 

0.04 

< 27 

8.81(0.01) a 

28(8) 

1100(360) 

> 1034 

76 

05358—mm3™ 

0.17(0.01) 

32(2) 

0.20(0.01) 

37(3) 

39.69(0.05) 

6.60(0.01) 

210(12) 

180(13) 

18 

G034-G2 

0.04 

< 6.95 

0.04 

< 6.8 

4.314(0.001)“ 

5.95(0.08) 

> 856 

> 872 

19 

G034-F1 

0.05 

< 6.94 

0.04 

< 5.8 

4.4(0.l) a 

3.9(0.3) 

> 566 

> 672 

11 

G034-F2 

0.04 

< 7.83 

0.04 

< 6.6 

9.20(0.03) 

1.53(0.01) 

> 195 

> 232 

18 

G028-C1 

0.04 

< 7.02 

0.04 

< 5.9 

5.735(0.006)“ 

8.5(0.1) 

> 1217 

> 1445 

15 

120293—WC 

0.063(0.006) 

17(2) 

0.08(0.01) 

20(3) 

7.150(0.003)“ 

11.0(0.8) 

65(69) 

550(98) 

29 

22134—G™ 

0.03 

< 5.6 

0.04 

< 6.6 

7.80(0.03) 

1.30(0.01) 

> 232 

> 197 

18 

22134—B 

0.04 

< 5.8 

0.04 

< 5.7 

2.83(0.07)“ 

1.8(0.3) 

> 316 

> 322 

14 

HMPO 

100117—MM 1 

0.05(0.01)* 

14(3) 

0.05 

< 12.8 

2.88(0.01)“ 

3.0(0.3) 

220(60) 

> 233 

27 

AFGL5142-MM 

0.072(0.009) 

32(4) 

0.042(0.003)* 

18(1) 

9.41(0.02)“ 

23.5(1) 

740(97) 

1300(210) 

53 

05358—mm 1 

0.17(0.02) 

61(6) 

0.18(0.02) 

65(8) 

36.63(0.05) 

11.78(0.02) 

190(20) 

180(23) 

43 

18089-1732 

0.08(0.09) 

22(9) 

0.10(0.01) 

27(3) 

9.312(0.002)“ 

21(2) 

1000(400) 

800(100) 

31 

18517+0437 

0.058(0.008) 

21(3) 

0.089(0.005) 

32(2) 

25.98(0.03) 

8.35(0.01) 

390(53) 

260(20) 

43 

G75—core 

0.06 

< 22 

0.05 

< 1.9 

13.95(0.06) 

4.49(0.02) 

> 209 

> 240 

43 

120293—MM1 

0.086(0.007) 

56(5) 

0.19(0.02) 

122(14) 

16.410(0.001)“ 

44.62(0.05) 

790(67) 

370(42) 

82 

121307 

0.05 

< 7.28 

0.05 

< 6.3 

2.13(0.05)“ 

2.2(0.3) 

> 301 

> 346 

10 

123385 

0.05 

< 19 

0.06 

< 23 

4.74(0.03) 

1.53(0.01) 

> 81 

> 66 

43 

UC Hu 

G5.89-0.39 

0.18(0.01) 

84(6) 

0.13(0.01) 

61(6) 

66.3(0.5) 

27.3(0.3) 

320(25) 

450(51) 

57 

119035—VLA1 

0.06 

< 31 

0.05 

< 24 

3.820(0.008)“ 

18(1) 

> 572 

> 754 

66 

19410+2336 

0.18(0.02) 

96(8) 

0.16(0.01) 

82(8) 

21.27(0.01)“ 

43(5) 

450(96) 

500(100) 

64 

ONI 

0.21(0.01) 

72(3) 

0.29(0.02) 

94(8) 

12.795(0.001)“ 

33(4) 

460(22) 

350(33) 

39 

122134—VL A1 

0.05 

< 24 

0.04 

< 18 

4.71(0.03) 

1.94(0.01) 

> 80 

> 108 

57 

23033+5951 

0.048(0.007) 

23(3) 

0.065(0.006) 

30(3) 

7.660(0.001)“ 

20(6) 

900(400) 

700(250) 

57 

NGC7538-IRS9 

0.05 

< 24 

0.05 

< 21 

14.93(0.04) 

6.13(0.02) 

> 255 

> 294 

57 


a lines with well-constrained optical depth, for which we do not list the integrated intensity 
but the parameter T ant x r t , from which we have derived the column density as explained in 
Sect. 13.11 

w ’’warm” HMSC; 
t tentative detection 

















Table 2: Integrated intensity of the C 15 N and 13 CN transitions, total beam-averaged column 
densities computed as explained in Sect. 13.21 and corresponding 14 N/ 15 N isotopic ratio. We 
do not list the upper limits of the integrated intensity of 13 CN and C 15 N(2-1) for the sources 
undetected in both transitions. To compute N( CN), the 13 CN column density was corrected 
by the 12 C/ 13 C ratio at the source galactocentric distance, D G q (Col. 5), assuming the 
relation between 12 C/ 13 C and D GG derived for CN by Milam et ah (120051) . Uncertainties 


have been derived as explained in the caption of Table [Tj for the optically thin lines. 


source 

C 15 

f T M b<1v 

K km s —1 

N(2- 1) 
V(C 15 N) 
(x 10 12 cm“ 

f T M b<1v 
~ 2 ) K km s -1 

13 CN(2 1) 

Dec N(OSS) 
kpc (xl0 14 cm -2 ) 

N (ON) 

JV(C 1B N) 

HMSC 

I00117-MM2 







AFGL5142-EC” 

0.14(0.02) 6 

1.6(0.2) 

0.35(0.02)“ 

10.3 

5.3(0.3) 

330(67) 

05358—mm3 1 " 

0.12(0.02) 6 

0.7(0.1) 

0.36(0.02)“ 

10.3 

2.7(0.1) 

400(90) 

G034-G2 







G034-F1 







G034-F2 







G028-C1 







I20293-WC 







22134-G 1 " 

0.06(0.02) 6 

0.3(0.1) 

0.16(0.02)“ 

9.5 

1.1(0.1) 

330(150) 

22134—B 







HMPO 

I00117-MM1 

< 0.06 

< 0.4 

0.13(0.02)“ 

9.5 

1.0(0.1) 

> 250 

AFGL5142-MM 

0.28(0.02)*’ 

2.4(0.2) 

0.40(0.02)“ 

10.3 

4.7(0.2) 

190(23) 

05358-mml 

0.12(0.02)*’ 

1.5(0.2) 

0.36(0.02)“ 

10.3 

3.7(0.2) 

240(40) 

18089-1732 

0.14(0.03)*’ 

0.9(0.2) 

0.80(0.02)“ 

5.1 

4.0(0.1) 

450(100) 

18517+0437 

0.10(0.02)*’ 

0.8(0.2) 

0.34(0.02)“ 

6.5 

2.4(0.2) 

310(80) 

G75-core 

0.19(0.03)*’ 

1.4(0.2) 

0.40(0.02)“ 

8.4 

3.4(0.2) 

240(50) 

I20293-MM1 







121307 

< 0.06 

< 0.4 

0.09(0.02)“ 

9.3 

0.7(0.2) 

> 175 

123385 







UC Hu 

G5.89-0.39 

0.14(0.03)“ 

3.5(0.8) 

0.12(0.03)“ 

7.2 

12(3) 

350(160) 

119035-VLA1 

< 0.06 

< 0.6 

0.16(0.02)“ 

7.0 

1.6(0.2) 

> 270 

19410+2336 

0.07(0.02)'* 

1.9(0.5) 

0.07(0.02)“ 

7.7 

8(2) 

430(250) 

ONI 

0.07(0.02)'* 

1.4(0.4) 

0.06(0.02)“ 

8.0 

5(2) 

350(220) 

122134-VLA1 

0.09(0.02)*’ 

0.8(0.2) 

0.18(0.02)“ 

9.5 

2.1(0.2) 

250(84) 

23033+5951 

< 0.07 

< 0.7 

0.21(0.02)“ 

10.3 

2.6(0.2) 

> 370 

NGC7538-IRS9 

0.08(0.02)*’ 

0.7(0.2) 

0.14(0.02)“ 

9.9 

1.7(0.2) 

230(90) 


“ ’’warm” HMSC; 

a from the hyperfme components J,Fi,F = (5/2, 3,4) — (3/2, 2, 3), (5/2, 3, 3) — (3/2, 2, 2) 
and (5/2, 3, 2) - (3/2, 2,1) ; 

b from the hyperfme components J, F — (5/2, 3) — (3/2, 2) and J,F = (5/2, 2) — (3/2,1); 
c from the hyperfme component J, F \, F = (3/2,1,1) — (1/2, 0,1); 
d from the hyperfme component J, F = (3/2, 2) — (1/2,1). 
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Fig. 1.— Left panels: spectra, from top to bottom, of 15 NNH + (1-0), N 15 Nff + (1-0) and 
N 2 H+ (1-0) towards 05358-mm3. The velocity interval shown is ±20 km s _1 from the 
systemic velocity. For the 15 NNH + line, the red curve represents the best Gaussian fit, while 
for the other lines corresponds to the best fit to the hyperfine structure. 

Right panels: spectra of C 15 N (N = 2 — 1, top) and 13 CN (iV = 2 — 1, bottom) towards 
05358-mm3 centred on the hyperfine components mentioned in Sect. 13.21 and indicated by 
the red vertical lines, the length of which is proportional to their relative intensities. 
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the NNH + / 15 NNH + and CN/C 15 N ratios can be measured in the same object, they have 
the same 14 N/ 15 N ratio, within the error bars, only in nearly half of the sources detected in 
both N 2 H + and CN isotopologues, as shown in Fig. [3J However, globally the lower range 
of variation of the CN/C 15 N could be a result biased by the fact that fewer sources were 
detected in C 15 N than in the N 2 H + isotopologues. Follow-up higher sensitivity observations 
will be helpful to address this aspect. 


Note that, since the CN/C 15 N has been derived by observations of 13 CN, there is also 
the possibility that an error is introduced because of the variation of the 12 C/ 13 C introduced 
by depletion of 13 C (see Sect. □])• However, the recent, very detailed study by Roueff et 
ah 1)20151 ) predicts variations up to ~ 20% (CN/ 13 CN ratio from 67 to 80), namely within 
the errors of our measurements. Thus, we conclude that our estimates of the CN/C 15 N are 
substantially correct. 


The measurements of both NNH + / 15 NNH + and CN/C 15 N are at odds with the recent 
model by Roueff et al. (120150 . which predicts that the 14 N/ 15 N should not vary in both 
species, unless the age of the condensation is shorter than 10 5 yrs. This possibility, however, 
seems to be ruled out by our observations, as the 14 N/ 15 N is not different in the three 
evolutionary groups. We conclude that some important routes of nitrogen fractionation are 
still missing in the models. 


A hint of what they could be might be obtained considering the D-fractionation of NNH + 
towards the same sources where we measured the N-fractionation. The bottom panels in 
Fig. [2] show NNH+/ 15 NNH+ and NNH+/N 15 NH+ as a function of the NNH+/NND+, as 
previously measured by Fontani et al. 12 0111) . Although not very strong, the data show a 
possible anti-correlation between the two parameters: the larger the 14 N/ 15 N, the lower the 
H/D. Non-parametric statistical tests confirm this trend: the Kendall’s r rank correlation 
coefficient is ~ —0.5 when the 14 N/ 15 N is derived from either 15 NNH + and N 15 NH + , and the 
Pearson’s correlation coefficient is ~ —0.6 and —0.5, respectively. Moreover, the probability 
for a chance correlation is below the 2% in both cases. Also, we have performed least square 
fits to the data, which provide slopes of —0.98 and —0.80 for the 14 N/ 15 N ratios derived from 
15 NNH + and N 15 NH + , respectively. Based on this suggestive anti-correlation, we conclude 
that the N-fractionation does not seem to be linked to the parameters that are known to affect 
the D-fractionation (if not to obtain the inverse effect): temperature and CO depletion (or 
density). In other words, the reactions leading to N-fractionation of NNH + and CN indeed do 
not have an energy barrier, as predicted by Roueff et al. (1201511 . Finally, and consistently with 
this last sentence, we do not find any correlation between the measured N-fractionation and 
the excitation temperature, contrary to the results of Wampfler et ah (2014), who claimed a 
tentative trend towards a temperature-dependent 14 N/ 15 N fractionation. However, because 
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the anti-correlation found is not strong, care needs to be taken in the interpretation of this 
result. 


Conclusions 


We have observed the J = 1-0 rotational transitions of 15 NNH + , N 15 NH + and N 2 H + , 
together with N = 2-1 transitions of 13 CN and C 15 N, towards 26 massive star forming cores 
in different evolutionary stages, where the deuteration fraction has already been measured 
(Fontani et al. 2011 b We find 14 N/ 15 N in CN between ~ 190 and ~ 450, and in N 2 H + 
from ~ 180 up to ~ 1300. It is the first time that high 15 N fractionations are found in 
N 2 H + . However, in the fractionation process, time does not seem to play a role because 
the 14 N/ 15 N ratio does not depend on the evolutionary stage of the source. We find a 
suggestive anticorrelation between D and 15 N fractions, consistent with findings in Solar 
System material, where N isotope ratios largely vary and only in a few cases they have 
been found to correlate with D/H ratios (Mandt et al. 2014). This lack of correlation is 
also consistent with the prediction of chemical models, which rule out any link between 
fractionation of N and the causes of D enrichment, namely temperature and density (or CO 
depletion). 
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